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ABSTRACT 

The temporal and spectral properties of ultraluminous X-ray sources (ULXs, Lx > 2 x 10 39 ergs 
s" 1 ) and bright X-ray sources (Lx > 3 x 10 38 ergs s _1 ) are examined and compared in two extremely 
different host environments: the old elliptical galaxy NGC 1399 and the young, starforming Antennae 
galaxies (NGC 4038/4039). ULXs in NGC 1399 show little variability on either long or short time 
scales. Only 1 of 8 ULXs and 10 of 63 bright sources in NGC 1399 are variable at a confidence level 
of 90%. On long timescales, the NGC 1399 sources are steadier than most Galactic black hole X-ray 
binaries, but similar to GRS 1915+105. The outburst duration of the NGC 1399 sources is about 20 
yrs, again, similar to that of GRS 1915+105. The bright X-ray sources in NGC 1399 may be black hole 
X-ray binaries with giant star companions similar to GRS 1915+105. The brightest ULX (PSX-1) in 
NGC 1399 is coincident with a globular cluster, shows a hard spectrum with a photon index around 
1.5, and has a nearly constant luminosity around 5 x 10 39 erg s . It may be an intermediate-mass 
black hole (IMBH) in a hard spectral state. In contrast to NGC 1399, the ULXs in the Antennae are 
all variable and a large fraction of the bright sources (9 of 15) are also variable. The variability and 
luminosity of ULXs in the Antennae suggest they are black hole high mass X-ray binaries accreting 
via Roche-lobe overflow. A flare with a duration of about 5 ks is found from Antennae X-42. The 
most luminous ULX, X-16, with a very hard spectrum (r = 1.0~1.3) and a luminosity which varies 
by a factor of 10, could be an IMBH candidate. 

Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X-rays: 
galaxies — galaxies: individual (NGC 1399, the Antennae (NGC 4038/4039)) 



1. INTRODUCTION 

Ultraluminous X-ray sources are point-like X-ray 
sources in external galaxies with positions off the nu- 
cleus and apparent isotropic luminosities in excess of 
the Eddington luminosity of a 10Mm black hole (~ 
2 x 10 ergs s ; for re view see iFabbiano fc Whitell20M 
iMiller fc Colberfl2004j) . ULXs a re often found in act ively 
starforming galaxies, e. g. M82 llKaaret et al.1 120011 and 
the Antennae galaxies l|Zezas etldTl2002al) . There are 
statistically more ULXs in starburst ga lax ies compared 
to spi ral galaxies l)Kilgard et alJ I2005T) . iKaaret et al.l 
( 2004) investigated three star burst galaxies and found 
X-ray sources were preferentially near star clusters but 
displaced with higher luminosity sources nearer to clus- 
ters. All these suggest that ULXs trace young stellar 
objects. 

Some surveys have reported hig her numbers of ULXs i n 
elliptical galaxies than in spirals l)Colbert fc Ptaki r2002'l. 
but the elliptical galaxies su ffer more contamination 
from bac kground sou rc es llLopez-Corredoira fc Gutierrez! 
120061) . Ilrwin et all (|2003H found that bright X-rav 
sources in elliptical galaxies have spectra consistent with 
Galactic X-ray binaries in the high state, but that the 
number of sources in the ULX range is consistent with the 
number of expected background sources. However, some 
elliptical galaxies contai n ULXs associated wi th globular 
clust ers, e.g. NGC 1399 fAngelini et afll200l and NGC 
720 l|,Ieltema et al.ll20?)3T) . which strongly suggests an as- 
sociation with the elliptical galaxy. 

NGC 1399 and the Antennae are cases, respectively, of 
an elliptical galaxy and a starforming galaxy that host 



a large number of ULXs. lAngelini et al.l ((2001) com- 
pared Chandra and Hubble Space Telescope (HST) images 
and found 70% of low mass X-ray binaries (LMXBs) in 
NGC 1399 were associated with globular clusters (GCs); 
the GC associated sources are brighter than field sources 
on average. The X-ray source population in NGC 1399 
may be unusual because a recent survey of LMXBs in 
six elliptical galaxies shows that, in general, there is 
no significant difference in l uminosity betwe en GC as- 
sociated and field LMXBs ilKim el alj |21)0:.i. Several 
ULXs in NGC 1399 are detected but their natures are 
still unclear; they co uld be IMBHs or multiple LMXBs 
l|Angelini et alJ l200l . X -ray sources in the Antennae 
have been well studi ed along with a series of observa- 
tions from Chandra (Fab biano et alJ 120011: iZezas et alJ 
2002a). Comparison with optical images shows that 
the Antennae X-ray sources are predominately associ- 
ated with young star clusters ijZezas et~aLll2002bt) . The 
Antennae ULXs have spectra consistent with ULXs in 
nearby galaxies and are not supernova remn ants because 
they a re variable and lack radio emission l|Zezas et alJ 
I2002b|) . IZezas fc Fabbiano! l)2002j) argued that IMBHs 
could account for a few of the ULXs in the Antennae 
but not most of them, which were more likely to be su- 
percritical accreting stellar mass X-ray bin aries due to 
their d isplacement from young star clusters (iZez as et alJ 
2002a). The ULX X-37 in the Antennae has be en iden- 
tified as a background quasar llClark et al.ll2005l) . 

The most interesting ULXs are potential IMBHs. NGC 
1399 and the Antennae are two extremely different host 
environments for ULXs. The comprehensive connection 
between ULXs and young stellar objects makes ULXs in 
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TABLE 1 

Chandra Observations of NGC 1399 and the Antennae. 



Galaxy 


index 


obs. ID 


exposure (ks) 


Instrument 


obs. date 




(a) 


320 


3.5 


ACIS-I/NONE 


1999 Oct 18 




(b) 


319 


55.9 


ACIS-S/NONE 


2000 Jan 18 




(c) 


239 


3.6 


ACIS-I/NONE 


2000 Jan 19 


NGC 1399 


(d) 


49898 


13.1 


ACIS-S/HETG 


2000 Juu 15 




(e) 


240 


43.5 


ACIS-S/HETG 


2000 Jun 16 




(f) 


2389 


14.7 


ACIS-S/HETG 


2001 May 08 




(g) 


4172 


44.5 


ACIS-I/NONE 


2003 May 26 



(i) 


315 


72.2 


ACIS-S/NONE 


1999 Dec 01 


(ii) 


3040 


69.0 


ACIS-S/NONE 


2001 Dec 29 


(iii) 


3043 


67.1 


ACIS-S/NONE 


2002 Apr 18 


The Antennae (iv) 


3042 


67.3 


ACIS-S/NONE 


2002 May 31 


(v) 


3044 


36.5 


ACIS-S/NONE 


2002 Jul 10 


(vi) 


3718 


34.7 


ACIS-S/NONE 


2002 Jul 13 


(vii) 


3041 


72.9 


ACIS-S/NONE 


2002 Nov 22 



NGC 1399 mysterious. A comparison of ULXs between 
these two system may shed light on their natures. We 
analyzed Chandra archival data for both galaxies (§ |2J 
in order to measure their long-term and short-term time 
variability (§ I2.1|) and spectral variation (§ 12. 2|) . We dis- 
cuss the results and the implications in § El We adopt the 
dista nce to NGC 1399 as 20.5 Mpc ijMerritt fc Ferraresd 
I2001D and the dista nce to the Antennae as 13.8 Mpc 
( baviane et alJ2004(l . This new distance to the Antennae 
is obtained by measurement of the red giant branch tip 
and appears accurate to within 2 Mpc, while the older 
distance was estimated from the recession velocity. The 
new distance is more likely to be accurate. We note that 
this revised distance to the Antennae has reduced the 
Antennae source luminosities by a factor of 4 relative to 
previous X-ray studies. 

2. CHANDRA OBSERVATIONS AND DATA REDUCTION 

Chandra has made seven observations of NGC 1399 
during 1999-2003 and seven observations of the Anten- 
nae during 1999-2002. A study of some o f the Antennae 
observations has previously been made (Fabbiano et al 
2003). The observations are listed in Tabled in which 
we assigned letters to index observations of NGC 1399 
and Roman numerals to observations of the Antennae. 
Only data from the Advanced CCD Imaging Spectrom- 
eter (ACIS) were used in this paper. Observations (a), 
(c), (g) were made with the ACIS-I and (b), (d), (e), (f) 
were from the ACIS-S but with the High Energy Trans- 
mission Grating Spectrometer (HETGS) applied on (d), 
(e) and (f). Observations from (i) to (vii) were all made 
with the ACIS-S detector. 

Level 2 events files were reduced from level 1 data us- 
ing CI AO 3.3 with CALDB 3.2.0. We created images 
and exposure maps for each observation and searched for 
point-like sources with the wavdetect command. Expo- 
sure maps were produced by assuming a source spec- 
trum of power-law form with a photon index of 1.5 
l|Angelini et al] 120011 IZezas et al.l 12002a) and a Galac- 
tic absorption along the line of sight of 1.34 x 10 20 cm 2 
for NGC 1399 and 3.83 x 10 20 cm 2 for the Antennae 
ijDickev fc Lockmanlll990j) . A power-law spectrum and 
a disk blackbody spectrum have similar luminosities in 
the 0.3-10 keV band, but will be quite different in the 



energy range up to 100 keV. Since our data are only 
effective below 10 kev, an approximation with a single 
power-law spectrum is adequate. In NGC 1399, we se- 
lected only sources that were located at least 20 arcsec 
off the nucleus and within the D25 circle (6.9 arcmin in 
diameter from the RC3 catalog). In the Antennae, we 
selected only sources in the active starforming region, 
which is bounded by a 1.4 arcmin radius circle. 

Since the absolute astrometry accuracy of Chandra is 
0.6 arcsec 1 , which is worse than the angular resolution, 
a correction for the aspect is necessary between observa- 
tions before the same region could be applied in different 
observations for source selection. For NGC 1399, we se- 
lected all detected sources (> 3ct) in the region of interest 
as defined above in different observations. Then we took 
observation (b) as a reference and searched in other ob- 
servations for matched sources that were located within 1 
arcsec of the selected source in (b). A pointing offset was 
calculated by averaging the offsets of all matched sources 
for each observation. A similar correction was applied to 
observations of the Antennae with observation (i) used 
as the reference. 

2.1. X-ray Source Time Variability 

To examine source variability between observations, 
luminosity-luminosity (L-L) diagrams are presented for 
both galaxies. First, we selected two observations 
and generated source lists using wavdetect. Then we 
matched sources and found their luminosities in both 
observations. For an unmatched source in the list, we 
placed the same elliptical region in the other observa- 
tion and calculated the net counts by subtracting the 
background level, estimated with wavdetect, from the 
total counts encircled by the ellipse. Then we derived 
the luminosity from the net counts normalized with the 
exposure map and the galaxy distance. 

There are three long observations for NGC 1399: (b), 
(e), and (g). The L-L diagram between observation (b) 
and (e) is presented on the top panel of Fig. ^ an( l that 
between (b) and (g) is on the bottom panel of Fig. ^ A 
lower threshold of 3 x 10 38 erg s^ 1 is placed to exclude 

1 radius of the 90% uncertainty circle, see 
http : //cxc . harvard . edu/cal/ASPECT/celmon/ 
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Fig. 1. — The 0.3-10 keV L-L diagrams for X-ray sources above 3 X 10 38 
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Fig. 2. — The 0.3-10 keV L-L diagrams for X-ray sources above 3 X 10 38 erg/s in the Antennae between observation (i) and all others. The 
dotted line indicates the lower threshold in luminosity, and the dashed lines enclose a region without variability considering a systematic 
error of s/2. 
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most neutron stars as well as faint sources. If the source 
is undetected with a significance of 2a, a 2a upper limit 
is shown. Fig. [21 shows L-L diagrams for the Antennae 
sources between observation (i) and all others. 

We emphasize that we calculate luminosity from count- 
ing rate assuming a fixed spectrum. Our quoted lumi- 
nosities are the counting rate multiplied by a roughly 
constant factor (the factor is not exactly constant be- 
cause different instrument modes were used for NGC 
1399 and the detector response changed slowly for the 
Antennae). Because we are investigating the variability 
of individual sources, our comparisons of luminosities is 
very similar to a comparison of counting rates. If the as- 
sumed spectra form for a given source is incorrect, then 
the luminosities will be incorrect, but all luminosities for 
any individual source would be shifted by the same fac- 
tor (if the spectral shape of an individual source varies 
between observations, then that would be evidence for 
variability). We choose to present luminosities instead 
of count rates, because the former are easier to compare 
with measurements of other sources. 

We note that an estimate of the true luminosity de- 
pends on the spectral form assumed and the energy band 
used. The bolometric luminosity of the sources is likely 
larger than our quoted values, since emission is likely 
present both above and below the 0.3-10 keV band. To 
determine the magnitude of the uncertainty on the lu- 
minosity in our chosen energy band, we fit each source 
that has enough photons with an absorbed power-law 
and measured its intrinsic luminosity (the photon index 
and column density are allowed to vary; see details in 
Section 12.2(1 . Then we fit the spectrum and obtained 
another estimate of the luminosity by fixing the absorp- 
tion to the Galactic value and the photon index to 1.5. 
The ratio between these two luminosities varies between 
0.6-1.2 for bright X-ray sources in NGC 1399 and 1-2 
in the Antennae. Therefore, a bilateral systematic error 
of y/2 is derived. We also checked with XSPEC that the 
systematic error of \/2 could cover spectral changes from 
r = 1.5 — 3 for a single power-law spectrum, kT — 1.5 — 3 
keV for a single disk blackbody spectrum, or a reason- 
able combination of them; the column density was fixed 
at the Galactic value for NGC 1399 but was varied from 
the Galactic value to 10 21 cm~ 2 for the Antennae. That 
means this systematic error is valid for black hole bina- 
ries at most spectral states that have been observed in 
our Galaxy. The systematic error is shown as two dashed 
lines in the L-L diagram. Sources lying between the two 
dashed lines are considered constant within the accuracy 
of our measurement. 

Lightcurves for the ULXs, which we define as having 
a luminosity in 0.3-10 keV band above 2 x 10 39 ergs s _1 
in at least one observation, and for the six most vari- 
able sources are shown in Fig. for NGC 1399 and in 
Fig. 0] for the Antennae. The variability is defined as 
the maximum ratio in luminosities, which could be a 
number with a standard error if the lowest luminosity 
is detected or a. 90%: lower bound if the lowest luminos- 
ity is an upper limit. The luminosities in the lightcurves 
are derived slightly differently than those in the L-L di- 
agram. First we generate a complete source list, which 
is based on detected sources in the reference observation 
((b) or (i)), as well as sources only detected in other 



observations. We place the same circular region for a 
given source in all observations and calculate the lumi- 
nosity. This difference in technique causes a mild dif- 
ference with the luminosities in the L-L diagrams but is 
consistent within the errors. There are 8 ULXs found in 
NGC 1399, PSX-1, PSX-2, PSX-4, PS X-7, PSX-3, PSX- 
21, P SX-25 and PSX-6 (names are from lHumphrev et alJ 
2003), sorted by the maximum luminosity. Among them, 
PSX-1, PSX-2, PSX-25 and PSX-6 are found to be as- 
sociated with GCs. The GC identifi cations are from the 
H ST image in lAngelini et ail (|2001|) and the source list 
of lDirsch et all l)2003D . There are 6 ULXs detected in 
the Antennae, X-16, X-ll, X-44, X -42, X-37 and X-29 
(names are from lZezas et al.ll2002a|) . sorted in the same 
way. We note that the Ant ennae X-37 has be en identi- 
fied as a background quasar llClark et al |2005l) and X-29 
is the nucleus of NGC 4039 l|Zezas et al.l l2002bL We ex- 
cluded these two ULXs as well as the nucleus of NGC 
4038 (X-25) from all analysis in our paper. 

For all the X-ray sources in NGC 1399 and the Anten- 
nae, we calculated their variability defined as described 
above. We defined the variability function for a popula- 
tion as the number of sources with a maximum variabil- 
ity larger than a certain value. The variability functions 
for all bright sources (Lx > 3 x 10 38 ergs s _1 in 0.3-10 
keV) are shown in the top panel of Fig. El The fraction 
of bright sources which are variable (with a confidence 
level of 90% above the systematic error of V2) is 10/63 
for NGC 1399 and 9/15 for the Antennae. The fraction 
of ULXs which are variable is 1/8 for NGC 1399 and 
4/4 for the Antennae. For NGC 1399, the only variable 
ULX is PSX-21. The color of the optical counterpart 
(C-R = 0.69, Washington C and Kron R. iDirsch etai] 
( 2003)) for this source is inconsistent with that of a glob- 
ular cluster, but the color may be affected by a defect in 
the i?-band image (T. Richtler private communication). 
Of the other 9 variable sources, only two have optical 
counterparts within 1". Therefore, the majority of the 
variable sources are not associated with globular clus- 
ters. This is in sharp contrast to the fact that 70% of all 
sour ces in NGC 1399 are associated with globular clus- 
ters (JAngelini et alJl200lf) . 

The detection threshold in all observations is well be- 
low our ULX luminosity threshold. However, the de- 
tection thresholds in the shortest observations of NGC 
1399 is close to our luminosity cutoff of 3 x 10 38 ergs s _1 
for bright sources. We examined two long observations 
spaced by similar periods, i.e., (b) and (g) for NGC 1399 
and (i) and (vii) for the Antennae, to check if the dif- 
ference in variability between NGC 1399 and the Anten- 
nae still persists when only observations with detection 
thresholds well below our luminosity cutoff are used. We 
found that 7 of the 63 bright sources in NGC 1399 and 
4 of the 15 bright sources in the Antennae are variable. 
This is consistent with the conclusions that bright X-ray 
sources in NGC 1399 are less variable than those in the 
Antennae. 

To compare with the variability function of X-ray bina- 
ries with known characteristics, we selected 14 Galactic 
black hole LMXBs 2 and 25 neutron star LMXBs 3 and 

2 GRO J0422+32, GRO J1655-40, GRS 1009-45, 
GS 1124-684, GS 2000+250, GS 2023+338, GX 339-4, V4641 Sgr, 
X0620-003, X1543-475, X1705-250, XTE J1118+480, 
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Fig. 3. — Lightcurves of all ULXs (top) and 6 most variable sources (bottom) in NGC 1399. The labels indicate the maximum change in 
luminosities, which is a number with a standard error or a 90% upper limit. The dates for observations (c) and (e) arc shifted to 30 days 
later for clarity. 
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(imax > 3 X 10 38 ergs s _1 ) in NGC 1399 and the Antennae 
(top). Their normalized variability functions are shown in the 
bottom panel to compare with variability functions of 14 black 
hole LMXBs, 25 neutron star LMXBs, the black hole LMXB GRS 
1915+105 and the neutron star LMXB, 4U 1820-30, in a Galactic 
GC. 



calculated their variability functions using data from the 
All Sky Monitor (ASM) onboard the Rossi X-ray Timing 
Explorer (RXTE). For each source, we got one variability 
value by sampling the ASM one day average lightcurve 
spaced in the same way as those seven Chandra obser- 
vations of NGC 1399. Then we shifted the first date 
to 10 days later and got another value until there were 
100 values. The variability of the source was derived 
from an average of these 100 values and then the vari- 
ability functions were obtained, respectively, for sources 
of black hole and neutron star LMXBs. We sampled the 
ASM data spaced in the same way as the Antennae ob- 
servations and the results looked similar. We also exam- 
ined the variability of two individual Galactic sources: 
GRS 1915+105 and 4U 1820-30. The Galactic black 
hole binary GRS 1915+105 has a luminosity in the ULX 
range and is unusual in t hat it has been a p ersistent X- 
ray emitter for 14 years ijDone et al.ll2004j) . It is also 
the most massive stellar mass black hole that has been 
found so far. The bright neutron star LMXB 4U 1820-30 
(=X1820-303) residing in the Galactic GC NGC 6624 
has a similar environment with sources in NGC 1399 

XTE J1550-564, XTE J1859+226 

3 Cir X-l, Cyg X-2, GS 1826-238, GX 13+1, GX 17+2, GX 3+1, 
GX 340+0, GX 349+2, GX 354-0, GX 5-1, GX 9+1, GX 9+9, 
LMC X-2, Sco X-l, Ser X-l, X0614+091, X1254-690, X1543-624, 
X1608-522, X1705-440, X1735-444, X1812-121, X1820-303, 
X1822-000, X2127+119 



and is an ultracompact binary, though its luminosity is 
low relative to the sources considered here. The vari- 
ability function for each of these sources was obtained 
from 100 ASM flux measurements. The normalized vari- 
ability functions of selected black hole and neutron star 
LMXBs, GRS 1915+105, 4U 1820-30, NGC 1399 and 
the Antennae are shown together in the bottom panel of 
Pig. 13 

To estimate the typical outburst duration of sources 
in NGC 1399, we used the two deepest observations of 
NGC 1399. These are (b) and (g) and are separated 
by AT — 3.35 yrs. If a source was detected in both 
observations, we assume the source lasted in its outburst 
state for AT. We found 60 bright sources and 8 ULXs 
had at least one detection in either (b) or (g) ,and among 
them 50 bright sources and 7 ULXs had two. Assuming 
a fixed outburst duration Tb for X-ray sources in NGC 
1399, we estimated that Tq — AT/(1 — rf), where rj was 
the fraction of twice detected sources among sources with 
at least one detection. We thus derived Tb = 20.1 yrs 
for bright sources and Tb = 26.8 yrs for ULXs. 

To search for short-term variability for every X-ray 
source in each observation, we used the K-S test to ex- 
amine if the 0.3-10 keV X-ray count rate is constant. 
No source in NGC 1399 shows any short-term variability 
above a significance level of 3<r. However, four Antennae 
X-ray sources, including two ULXs, are variable with a 
significance level above 3er. A 3.8ct probability of non- 
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Fig. 6. — Lightcurves for the Antennae X-16 in the observation 
(v) and X-42 in the observation (vii) and the PSD for X-42. The 
top panel shows lightcurves in 1-7 keV band binned with 20 ks. The 
bottom panel presents the PSD in the 1-7 keV band calculated from 
a 512-point FFT logarithmically binned by a factor of 1.2. The 
dotted line indicates the Poisson level, and points below which are 
plotted with 2a upper limits. 
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TABLE 2 

Spectral parameters with 90% errors of PSX-2 and 
PSX-4 in NGC 1399 



source 


model 


n H (cm 2 ) 


kT (keV) or V 


X 2 /dof 


PSX-2 


diskbb 
powerlaw 


l-3±g"5 x 10 20 
6.1±\i x 10 20 


40+ - 02 
9 7+O.2 
z - '-0.1 


30.6/29 
37.7/29 


PSX-4 


diskbb 
powerlaw 


1.3+11 x 1020 
i.2+il x 10 20 


42+ - 03 

9 q+0.3 
z '°-0.2 


18.7/26 
20.2/26 



constant emission is shown for the Antennae X-16 in ob- 
servation (v), and an 8.0a level for X-42 in observation 
(vii). The total number of trials, equal to the number 
of detected sources multiplied by the number of observa- 
tions is 105. Taking into account the number of trials, 
the variability from X-16 is significant at the 98.7% level, 
while the variability from X-42 is highly significant with 
a chance probability of occurrence of4.1xlCP 7 . We ex- 
amined the light curves in several bands and found that 
the variability is most pronounced in the 1-7 keV band; 
there is no significant variability below 1 keV. The 1-7 
keV lightcurves for X-16 and X-42 and the power spec- 
trum density (PSD) for X-42 are shown in Fig. El Strong 
variability at a time scale around 20 ks in X-42 is explicit. 
The PSD looks like a flat form that breaks at 4 x 10~ 4 
Hz, which is similar to the break PSDs of Galactic black 
holes. However, the power spectrum is dominated by a 
single flare with a duration of about 5 ks about halfway 
through the observation. Therefore, it would be impru- 
dent to over interpret the power spectrum and attempt 
to derive a mass for the compact object. We note that 
no significant power is detected in any other observation. 

2.2. X-ray Spectra for ULXs 

We examined all the X-ray spectra for ULXs in both 
galaxies with XSPEC 11.3.1. All sources in NGC 1399 
and the Antennae could be reasonably described by an 
absorbed power-law, whereas some sources in the Anten- 
nae may show line f eatures but with low significance (see 
iZezas et all l|2002aj) ). Two sources (PSX-2 and PSX-4) 
in NGC 1399 present an alternative favored fitting with 
a multicolor disk (MCD) blackbody (diskbb in XSPEC) 
model (see Fig. |7J| . All fittings adopt an absorption with 
solar abundance (wabs in XSPEC) as a free parameter 
but setting the lower bound to the Galactic value. 

Fitting parameters of sources PSX-2 and PSX-4 are 
shown in Table The best fitted 0.3-10 keV luminos- 
ity versus photon index in every observation that con- 
tains enough photons are plotted in Fig.|8]for NGC 1399 
ULXs and Fig. H for the Antennae ULXs. The best fitted 
column densities for ULXs in NGC 1399 are consistent 
within 3(J with the Galactic value, except for PSX-2 in 
observation (b) which differs by nearly 6c The best fit- 
ted column densities for ULXs in the Antennae are shown 
on each plot panel in Fig. ED All photon indices and lu- 
minosities in Fig. |S| & are obtained by fitting a single 
power-law model with absorption. 

In Fig. |H1 a weak correlation between the luminosity 
and the power-law photon index is found for all ULXs 
except PSX-1 in NGC 1399 with a correlation coefficient 




PSX-2 



0.5 1 2 

channel energy (keV) 




channel energy (keV) 

Fig. 7. — X-ray spectra of the source PSX-2 (top) and PSX-4 
(bottom) in NGC 1399 from observation (b) fitted by a diskbb 
spectrum. The referred disk temperatures are 0.40 keV and 0.42 
keV, respectively for PSX-2 and PSX-4. However, both data sets 
can also be fitted by a power-law model with a higher absorption. 
See Table |21 for all fitting parameters. 



of 0.66 and a probability of 3.7% for a coincidence from 
uncorrelated data. This correlation does not exist in all 
ULXs in the Antennae, but in observations for individual 
sources X-ll, X-44 and X-42, the correlation coefficients 
with chance probabilities in brackets are 0.61(14.6%), 
-0.57(17.5%) and -0.66(10.6%), respectively. These 
correlations are weak and more observations are needed 
to expand the samples. 

3. DISCUSSION 

The most dramatic difference between NGC 1399 and 
the Antennae is that the ULXs as well as the bright X- 
ray sources in the elliptical galaxy are less variable. In 
Fig. m only a couple sources show apparent variability 
in the L-L diagram of NGC 1399 between (b)-(e) or (b)- 
(g) observations. In Fig. |3 the number of sources in the 
Antennae with significant variability is obviously larger 
than that in NGC 1399. This is verified by lightcurves 
of ULXs and the variability functions of bright sources. 

The variability functions of Galactic LMXBs show that 
most black holes are highly variable, with the notable ex- 
ception of GRS 1915+105, whereas bright neutron stars 
are less variable. The normalized variability function 
(bottom panel of Fig. 03) presents the bright sources in 
NGC 1399 as similar to 4U 1820-30, a little less vari- 
able than the group of neutron star LMXBs and GRS 
1915+105, and far less variable than the group of black 
hole LMXBs. ULXs in NGC 1399 are even less variable 
than other bright sources, with 7 of 8 showing constant 
flux and the other having little variability. None of the 
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bright X-ray sources or ULXs in NGC 1399 show short- 
term variability with a significance level of 3tr in any indi- 
vidual observation. The lack of variability should provide 
a strong constraint on the nature of these sources. 
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Fig. 8. — 0.3-10 keV luminosity vs. power-law photon index (T) 
for the 8 ULXs in NGC 1399. All points in the plot come from 
observations that have enough photons for a spectral fitting. PSX- 
1, PSX-2, PSX-4 and PSX-7 have two points from observations (b) 
and (g), and rest sources have one from the observation (b). The 
absorption in the fitting is set as a free parameter but with the 



lower bound to the Galactic value (1.34 X 10 



1 ). The best 



fitted nn are all close or equal to the Galactic value (less than 3cr 
except PSX-2 in observation (b), which is nearly 6cr). 

The long term variability also allow us to estimate the 
outburst duration of X-ray sources in NGC 1399. We 
found that the outburst duration is around 20 yrs for 
both bright X-ray sources and ULXs. 

A natural way to produce low variability would be if 
the bright X-ray sources in NGC 1399 consist of multiple 
unresolved LMXBs. However, a large number of sources 
is required. From the ASM data, we estimated that neu- 
tron star LMXBs are above 10% of their maximum flux 
80% of the time, while black holes are above their 10% 
of their maximum flux only 15% of time (the estimate 
for black holes is an upper bound because several of our 
selected black hole systems have had only one outburst 
during the lifetime of RXTE). Using the average lumi- 
nosity of 4U 1820-30 of - 5 x 10 37 ergs" 1 , about 8 neu- 
tron star binaries would be needed to produce the least 
luminous of our bright sources and 50 would be needed 
to produce the least luminous ULX. Large numbers of 
neutron stars are known to be present in some globular 
clusters in our galaxy, e. g. there are around 50 in the 
Galactic GC NGC 6266 iTPoolev et a,l1l200l . but only 
very rarely does a cluster contain more than one bright 
X-ray source. It is very puzzling to understand how in- 
dividual clusters in NGC 1399 could contain so many 
bright X-ray sources. Though the luminosity of black 
holes in the accreting phase is, generally, several times 
larger than that of neutron stars, the smaller duty cycle 
partially offsets the increased luminosity and comparable 
numbers of black hole binaries are needed to reach the 
same luminosity. 

Perhaps, the most promising interpretation of the 
NGC 1399 bright sources is that they are are individual 
black holes of a type similar to GRS 1915+105. As noted 
above, GRS 1915+105 has been a pe rsistent X-ray source 
since its discovery 14 years ago (|Castro-Tirado et al.l 
Il994|) and has reached luminosities in the ULX range. 
The variability function of GRS 1915+105 is similar to 
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Fig. 9. — 0.3-10 keV luminosity vs. power-law photon index 
(r) for the 4 ULXs in the Antennae in seven observations. The 
absorption in the fitting is set as a free parameter but with the 
lower bound to the Galactic value (3.83 X 10 20 cm -2 ). The best 
fitted nn values are consistent in all observations and the weighted 
average number is shown for each source. 



that of the NGC 1399 bright sources. The companion 
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star i n GRS 1915+105 is a K or M giant ijGreiner et al.1 
l2001|) . The time scale for the evolution of such a binary 
system to the mass-transf erring; phase is roughly 6 Gyr 
(Belczvnski & Buh]J l2002|) . Allowing for different com- 
panion star masses, similar mass transferring binaries 
could currently be present in either the old (~ 11 Gyr) 
or young (~ 2 Gy r) globular c lusters in NGC 1399 
l)Forbes et al J 12001^ . King et all l)2002j) has previously 
suggested that t he ULXs in ellipt i cals m ay be similar to 
GRS 1915+105. iPiro fc Bildster] pjp^ suggested that 
the bright X-ray sources in ellipticals are X-ray bina- 
ries with Roche-lobe filling giant companions with mass 
transfer driven by nuclear evolution of the companion. 
They predicted outburst durations similar to what we 
find for NGC 1399. For the sources associated with glob- 
ular clusters, another possibility is that the companion 
stars are initially isolated stars which are captured by the 
black hole after evolving to the giant phase. The capture 
cross-section is enhanced for stars in the giant phase. 

However, GRS 1915+105 shows strong variability on 
short time scales, which we do not observe from the 
bright X-ray sources in NGC 1399. The sensitivity of 
Chandra to short timescale variability is limited due 
to the low counting rate, but our detection of short 
timescale variability from sources in the Antennae (with 
comparable counting rates) indicates that similar flares 
would have been detected from the NGC 1399 sources. 
The short timescale variability of sources in NGC 1399 
should be investigated in more detail with new observa- 
tions; 

iBildsten fc Delovel l)2004|) has suggested that the bright 
X-ray sources in NGC 1399 may be neutron stars 
accreting from white dwarf companions in ultracom- 
pact binaries, similar to 4U 1820—30. The variabil- 
ity of 4U 1820-30 is similar to the NGC 1399 bright 
sources. Hence, neutron star ultracompact binaries 
may account for some or all of the bright sources if 
mildly super-Eddington (apparent) luminosities can be 
achieved. However, the short lifetimes expected for the 
high accretion rates required given the very low-mass 
companions (0.04 — 0.08M Q ) make it unlikely that the 
ULXs in NGC 1399 are neutron star ultracompact bi- 
naries. Another possible interpretation of the bright 
sources is that they are neutron-star binaries similar to 
the neutron-star X-ray binaries known as Z sources. Z 
sources are persistent with variability similar to that 
seen from the NGC 1399 sources. Again, mildly super- 
Eddington (apparent) luminosities would be required to 
produce the dimmest of the bright sources and such ob- 
jects are unlikely to explain the ULXs . Interestingly, 
iWebbink. Rappaport. fc Savoniiel 1)1983!) suggested that 
the Z sources have giant star companions with mass loss 
driven by nuclear evolution. This would make the sys- 
tems the neutron star analogs of black hole/giant star 
binaries such as GRS 1915+105. 

The three most luminous X-ray sources in NGC 
1399 are interesting. The brightest one, PSX-1, has 
a hard spectrum (r = 1.2 ~ 1.5), little variability 
(imax/imin = 2.2 +_ 0.8), and a high luminosity around 
5 x 10 39 ergs s _1 . The spectrum indicates the source 
is not in its high/soft state (r = 2.1 — 4.8) but at the 
hard extreme of the low/hard state (r = 1.5 — 2.1) 
ijMcClintock fc Remillardll2006j) . The low variability in- 
dicates that the source is probably not relativistically 



beamed ijKording et alJ l2002|) . because relativistically 
beamed sources are expected to be highly variable. This 
source is associ ated with a bright G C with C — R = 2.12, 
rriR = 21.02 (Dirsch et alJ l2003j) . which means it is 
unlikely to be a background AGN. A huge number of 
sources would be required to explain the emission as a 
collection of neutron star or black hole X-ray binaries. 
If the source is an IMBH in the low state, its practi- 
cally constant luminosity is reasonable if compared to 
the long-term behavior of M 82 X-l, which varies i n lu- 
minosity by a factor of 2-3 ijKaaret et alJl2006albf ) and 
shows a hard spectru m with a power-law index r ps 1.6 
lIKaaret et al jl2006hT) . 

NGC 1399 PSX-2 and PSX-4 display typical behaviors 
of black holes in the high/soft state. PSX-2 is associ- 
ated with a bright GC (C - R = 2.24, m R = 21.36) 
but PSX-4 is not. The plausible cool disk (kT ps 0.4 
keV) emission is consistent with the soft spectral index 
(r ps 2.5). Their low variability (L max /L min = 2.5 ± 0.7 
for PSX-2 and 1.5 ± 0.6 for PSX-4) is suggestive of some 
Galactic black holes in the high state, but no Galac- 
tic source has been found to remain in the high state 
for four years <|McClintock fc Remillardl l2006'). The cool 
disk h as been a useful probe in finding IMBH candi- 
dates JColbert fc Mushotzkvl 119991 iKaaret et all l2003t 
iMiller et alJl2004HFeng fc Kaaretll2005D . However, their 
disk temperatures are somewhat higher than most IMBH 
candidates, and their luminosities (2 ~ 3 x 10 39 ergs 
s" 1 ) are also around the lower bound of t hose candidates 
({Miller et alJl2004HFeng fc Kaaretll2005D . These sources 
do not have high mass companion stars to p roduce geo- 
metric ally beamed emission as suggested bv lKing et alJ 
(2001). Therefore, it could be possible that these two 
sources are 10^100M Q black holes emitting at the high 
state if the disk components are true. One possible way 
to determine the reality of the disk emission is an accu- 
rate measurement of the column density, since these two 
sources, in particular PSX-2, show significantly higher in- 
trinsic absorption if modeled by a single power-law spec- 
trum. 

The bright X-ray sources in the Antennae are vari- 
able. This suggests that they are black hole binaries and 
accreting from massi ve stars to produce unstable disks 
ijKalogera et alJ2004j) . They are not as variable as black 
hole LMXBs (Fig. 0). Considering the variability to- 
gether with their high luminosity, we suggest they are 
HMXBs accreting via Roche-lobe overflow. There is no 
disk emission found in any of the 4 ULXs. 

The most luminous source, X-16, manifests a very hard 
spectra (T = 1.0 — 1.3) and has strong variability with a 
change of luminosity of ~10 (Fig. EJ). The source could 
be an IMBH with a mass of a few hundred to thousand 
Mq emitting at the hard state. The X-ray behavior is 
also consis tent with a relativi stically beamed jet pointing 
at us, but Ze zas et "all l)2002tJl pointed out that there is 
no radio emission found from any ULX in the Antennae. 
It is interesting that none of the ULXs shows the very 
soft spectrum normal ly seen in ULXs in nearby galaxies 
ijFeng fc Kaaretll2005|) . We note that X-37 was reported 
to have a soft component in its spectrum, but has been 
subsequently identified as a background AGN. X-ll and 
X-42 are likely to be in the hard state (r < 2) , X-44 looks 
like it is between the hard and soft states (r ~ 2). A 
plausible correlation (r = 0.61) is found in X-ll between 
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luminosity and T, which shows similar behaviors as stel- 
lar mass black holes, but we should be cautious because 
of the weak correlation coefficient. The X-ray properties 
of these three sources are consistent with HMXBs, but 
whether they are IMBHs or ~20M Q black holes emitting 
in near or super Eddington luminosities is still unclear. 
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